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Abstract 

We present FDTD calculations for transmission of light and other 
electromagnetic waves through periodic arrays of slits in a metallic 
slab. The results show resonant, frequency dependent, transmittance 
peaks for subwavelength widths of the slits which can be up to a 
factor of ten with respect to those out of resonance. Although our 
conclusions agree with previous work by Lezec and Thio as regards 
both the magnitude of the enhancement and the lack of contribution 
of surface plasmon polaritons of the metal surface to this effect, we 
derive an interpretation from a theory that deals with emerging beam- 
Rayleigh anomalies of the grating, and with Fabry-Perot resonances 
of the perforated slab considered as an effective medium. 
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1 Introduction 



Extraordinary optical transmission through hole arrays of subwavelength di- 
ameter in metal films has been subjected to extensive study since it was first 
reported the detection of transmitted intensity large enhancements, which 
were polarization dependent [1] and attributed to excitation of surface plas- 
mon polaritons (SPPs) of the film surface [2], [3]. Nevertheless, recently 
both the enhancements of the transmission and the role of the SPPs in the 
phenomenon have been reconsidered [4], [5] showing that neither the en- 
hancements are as large as initially proposed, nor the SPPs play a key role 
since the same observations are obtained in films of materials like Cr which 
has a very large permittivity imaginary part thus defining a very dissipative 
and broad plasmon. The same occurs in W and Si films that have a positive 
real part of the dielectric permittivity in the visible, and which of course do 
not support SPPs. These experiments arose a certain controversy as regards 
the understanding of the underlying physics behind phenomenon, and have 
thus warranted a closer look at the interpretation of these transmittance en- 
hancements as in a recent publication [6], where it is shown that the initially 
reported enhancements of several orders of magnitude are so when the ra- 
tio is taken with respect to the transmitted intensity obtained from Bethe's 
theory [7]. However, Bethe's approximation, only valid for apertures much 
smaller than the wavelength and that thus assumes the field as constant in- 
side them, is not applicable to this configuration and thus the aforementioned 
ratios are misleading. More realistic is to take the ratio of the resonantly en- 
hanced intensities with respect to those obtained, both in experiments and 
exact calculations, out of resonance as done in [4] and [6]. Then, these en- 
hancements are larger than without resonant transmission not by orders of 
magnitude but by a factor that in data obtained so far is not larger than 10. 

In this paper we present an interpretation of the enhanced transmission 
of light through a periodic array of slits in a metal slab by means of a model 
that combine an approximate analytical calculation of scattering from cor- 
rugated metallic interfaces with that of transmission through metal slabs. 
The results agree qualitatively with those of [4] and [6] for circular holes. 
We conclude that transmission peaks are produced either by those kind of 
Wood anomalies known as Rayleigh anomalies (RA), namely those due to 
the emerging beams appearing in the passing-off of a new diffracted order, 
[8], [9], [10], [11], and that are well known to giving rise to strong discontinu- 
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ities of the diffracted intensities, or as a result of the Fabry-Perot (FP) type 
resonance in the film layer as if it were homogeneous with an effective refrac- 
tive index. The FP resonance comes from reflection and transmission at the 
entrance and the exit interfaces of the slab. We shall show results both by 
an exact FDTD calculation and by means of a theoretical model that shows 
the physical meaning and nature of the enhancement phenomenon in the 
slit grating transmittance. In this connection, it should be pointed out that 
there is agreement that this kind of resonances appear in these slit arrays as 
recently shown in [12], [13], [14]; also related are those so-called horizontal 
and vertical surface plasmon resonances by the authors of ref.[15]. Other re- 
searchers like to interpret them in terms of coupling of SPPs at both surfaces 
of the film and of cavity modes in the slits [16]. Further studies along these 
lines can also be found in [17]. In addition, the interpretation of waveguide 
mode resonances and diffraction as responsible for the extraordinary trans- 
mission is pointed out in [18] and [19] where it is observed that transmittance 
of subwavelength metallic gratings may be nearly zero for frequencies corre- 
sponding to SPP excitation. Also, the introduction of film modulations to 
obtain similar resonant effects can be found in [20]. 

We shall first show how one slit does not transmit for s polarization but 
it does so for p-polarization. Then we shall argue that the transmission 
enhancement is not due to the SPP of the metal surface because with the 
slits present in the film this polariton is too extended and thus only one 
of its many Fourier components matches thus being not enough to yield the 
enhancement . By contrast, the Rayleigh anomaly is by contrast very sharply 
defined: that wavelength at which the z-component of the diffracted beam 
wavevector vanishes. 

As for the FDTD calculations [21], we address grids with a space step 
5x = AOnm and a time step 5t = 0.15x/c, where c is the light velocity. This 
sampling yields good convergent solutions. Periodic boundary conditions 
where used along the boundary parallel to the main direction of light prop- 
agation. On the other hand, the uniaxial perfectly matched layer (UPML) 
absorbing boundary condition, was used on the boundary perpendicular to 
light propagation; (showing no significantly different results for the particular 
structure here used to previous calculations done by means of Liao's condi- 
tions). Concerning the materials in these calculations, and as far as Ag is 
concerned, where SPPs have long been studied both in flat and periodically 
corrugated surfaces [22], [23], [24], and for which the imaginary part e 2 of its 
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permittivity e = e\ + is small enough so that the SPP is well defined [22] 
in Drude's model with a bulk plasmon wavelength \ p = 325nm u p ~ 3.8eV, 
it is known that for example in a periodic Ag interface with a small sin- 
gle Fourier component, the enhancement of the field due to the SPP may 
reach large values (100 times that of the incident field) [23]. However, when 
the amplitude of this harmonic increases, or many Fourier component exist, 
(this happens in e.g. a random rough interface [11] or in a step-like profile of 
the surface like in the slits addressed here), the enhancement is drastically 
reduced due to the enlargement of the SPP linewidth. In our calculations, 
the metal is characterized in Drude's model e(a>)/e/ = 1 — {uo 2 /uo 2 + 2iu5)) 
by fitting its permittivity e/, the bulk plasma frequency uo p and the damping 
constant 5 to the experimental data by Johnson and Christy [25] for Ag. The 
chosen parameters are tf = 6.8, uo p ~ 3.8eV and 5 ~ —0.02eV. 

2 Calculations 

We consider a plane electromagnetic wave normally incident on a lamellar 
grating consisting of 1-D slits as indicated in Fig.l, which we shall charac- 
terize by a profile function z = D(x) with respect to the entrance plane of 
the film: z — 0. This configuration was studied in [2]. The parameters are: 
period d = 3.5/im, slit width a = 0.5/xm; the thickness of the film h varies 
from 0.2/im to 4fim. 

We shall write the intensity transmitted by the grating in the far zone as: 

I(Q m ) = -\A(Q m )\ 2 , (1) 
% 

where the incident and diffracted wavevectos are: ko = (Ko,qo), k m = 
( Qm,Qm), Ki = k o = ( 27r A) 2 - Q m = K + 2irm/d, m = 0,1,2,..., q m = 

\Ao — Qm- 

The diffracted amplitude transmitted by the grating is written as [26]: 

A(Qm) = — f d dx'F{x') exp[-t(Q m x f + q m D{x'))\. (2) 
q m Jo 

In Eq.(2) the source function F(x) is given by the boundary conditions at 
the grating surface. On the other hand, for normal incidence: Q m = 27rm/d. 
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In order to explicitely compare the contribution of the FP resonance as 
compared with that due to the RA, we shall write F(x) as: 

F(x) = exp[iq D(x)]Vf exp(i^), (3) 
(j) being a certain phase and T being the film trasmittance: 



(1 - r) 2 1 + 4r/((l - r) 2 sin 2 (5/2)) ' V ; 

where 8 = q n e ffh + 2ip, n e ff is an effective index of the film averaged 
from the opaque and transparent portions of the slab, which makes sense for 
small wavelengths, and ip denotes a phase shift between the incoming and 
the reflected wave. 

In this way, the diffracted intensity for a given Q m becomes: 

I(Q m )=T\B(Q m )\ 2 . (5) 

B(Q m ) is obtained from Eqs.(2) and (3). We then perform the integration 
for the profile D(x) corresponding to the slit grating. The result is: 

B{Qm) = 4« — sm(— — ) sm[ J. (6) 

Qm & ^ 

The transmittance is next calculated by integrating I(Q m ) over all diffrac- 
tion orders, namely, at all angles of transmission. Figs.2(a)-2(d) correspond 
to the case dealt with in [2] . Both FDTD calculations (solid line) and results 
from Eq.(5) (broken line) are plotted. The slits are considered as practiced 
in an Ag slab of thickness h which varies from 0.2/xm to 4/im, the period 
being d = 3.5/im. Of course, the agreement between calculations and the 
analytical theory of Eqs. (2)-(6) is qualitative: these equations help us to 
interpret the origin of the resonant peaks of the transmission intensity, even 
though the numerical value of their results does not coincide with the exact 
one provided by the FDTD calculation. 

We observe peaks under p-polarization, (namely with the H vector along 
the y-axis in Fig.l). In this case there is no cut-off for transmission because 
the first mode of the aperture cavity has the index (1,0); namely, it corre- 
sponds to a homogeneous mode. This is the underlying reason why there 
is larger transmission for p-polartization. Conversely, under s-polarization, 
there is a frequency cut-off and thus an evanescent transmitted mode through 
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the aperture cavity which has the index (1,1), and hence there is no trans- 
mitted intensity since this cut-off exponentially dampens its intensity. (An 
illuminating discussion of this effect is presented in Jackson's text [27], Sec- 
tion 8.4). Only by overcoming this cut-off at high frequencies a transmitted 
intensity arises in s-polarization. 

On the other hand, in these figures, which show agreement with the cal- 
culation of [2], the two peaks, (in particular in Fig. 2(a) corresponding to 
h = 3/im and n e ff = 1.2), are well separated from each other. Neither 
of these peaks arise from surface plasmon polaritons (SPPs). The broader 
peak in the right is a FP resonance due to the contribution of T in Eq.(5). 
The narrower peak, placed at shorter wavelengths in Fig. 2(a), is a Rayleigh 
anomaly, namely one due to the passing-off of a new diffracted mode; as such 
it occurs at a sharply well defined value of the wavelength, namely that at 
which q m = 0. As the slab becomes thinner, the FP resonance peak shifts to 
smaller wavelengths and tends to overlap with the Rayleigh-emerging beam 
resonance. This is shown in Figs. 2(b)-2(d) corresponding to h = 2/j,m, 
h = 1.2/im and h = 0.6/im. In fact, Fig.2(d) shows that in the latter case 
the peak resonance is mainly from the pure emerging beam. In the cases of 
Figs. 2b) and 2(c), there is an overlapping of the FP and the emerging beam 
resonances with a subsequent splitting of both peaks. Therefore, this combi- 
nation of the FP plus the emerging beam resonances explain all the structure 
of the frequency dependent transmittance curves given by the exact FDTD 
calculations, in numerical agreement with those of [2]. 

It should be pointed out that, although the resonance frequency of SPPs 
of the flat Ag slab is very close to that of the RA of this grating, the large 
corrugation introduced on its surface by the presence of the slits, in addi- 
tion to to its broad Fourier spectrum as pointed out before, broadens the 
scattered field angular spectrum which, consequently, only weakly interacts 
with the SPP of the film flat surface. This is well known from the theory 
of surface-relief grating diffraction [23], [11], [24], where it is shown how the 
absorption of light by a grating, due to the SPP excitation, decreases as its 
depth increases. In addition, as we shall see later, these RA peaks equally 
appear for a non metal slab with slits as Cr. 

As the slab thickness h increases, more FP resonances appear at longer 
wavelengths. This is shown in Fig. 3 (a) for h = 4/im which corresponds to 
the calculation using Eq.(5) and that now exhibits two of these FP peaks: 
one near 5000nm and an additional broader one at lOOOOnm. As remarked 
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above, none of these peak structures appears in s-polarization; this is seen in 
Fig. 3(b) by the FDTD results. In this last figure we show the transmission 
curve slowly rising from the wavelength of 1.5/im and more abruptly towards 
the left to lfxm, which is precisely double of the slit width a = O.bfim, and 
thus constitutes the cut-off wavelength. This further demonstrates that the 
observation of transmission in these structures under p-polarization only, is 
not caused by SPP resonances, but is rather due to the above mentioned 
transmission characteristics of the waveguide that constitutes the aperture. 
We would also like to stress that these transmission enhancements, obtained 
at resonance, are not orders of magnitude larger, but given by factors not 
larger than 10. 

The main contribution of both FP and Rayleigh anomalies to this trans- 
mission behavior, and not of SPPs, is further shown in Fig. 3(c) which cor- 
responds to a slit grating perforated in a Cr slab of thickness h = 0.6/xm. 
Namely, identical to that dealt with in Fig. 2(d) except for the material that 
now has a completely different dielectric constant. In fact, its imaginary 
part that defines the sharpness of the SPP at the frequencies of Fig. 3(c) is 
10 times larger than for Ag. The transmitted intensity exhibits not only 
the same RA resonance, but also an almost identical transmission curve as 
that of the Ag grating at all frequencies, as one sees on comparison of Figs. 
2(d) and Fig.3(c). This result agrees with previous reports stating that these 
transmission peaks are equally obtained in slits practiced in metals and in 
dielectrics [4], [5], [28], [29]. 

On the other hand, as mentioned in Section 1, SPPs are well defined in 
Ag and Au in the visible . Then the width of the plasmon is small and its 
lineshape is sharp, (see [22]) . However, at the frequencies at hand the values 
of the imaginary part of Ag and Cr are about 5-35 and 44-200, respectively. 
This produces, (even more for Cr), a quite broad, and hence badly defined, 
SPP, if any. In addition, a SPP is a well defined harmonic wave. This is so 
when the roughness of the suface, (i.e. the h parameter in our grating) that 
holds the plasmon, is low. This leads to large enhancements in the near field 
intensity. However as the roughness increases, the plasmon width increases, 
and then the field enhancement is drastically reduced [23], [11]. Hence, SPPs 
cannot significantly contribute to these transmission peaks. On the contrary, 
for the Rayleigh resonances obtained in these calculations, the trasmitted 
intensity increases with the roughness h up to a certain value, like when the 
corrugation is a harmonic function [10], [11]. 
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3 Conclusions 



We have put forward an analytical approximate model, validated by FDTD 
calculations, that accounts for the transmitted intensity enhancements of 
light and other electromagnetic waves transmitted through a grating con- 
sisting of slits in a slab. The model has an interpretative value as it shows 
that these enhancements are due to both Rayleigh anomalies of the grating 
and to Fabry-Perot resonances of the slab. The contribution of surface plas- 
mon polaritons to this effect is negligible, although when the apertures have 
other geometry, like circular holes, their morphological resonances (localized 
plasmons, not to confuse with SPPs) do also play an important role in these 
enhancements [6] . On the other hand, the enhancements obtained are a fac- 
tor not larger than ten, but not by orders of magnitude, larger than out of 
resonance. We notice that there are critical experimental data that observed 
similar enhancement factors for Ag, W and Si in gratings of holes [4] , [5] in 
the visible, where the last two materials are dielectric, (namely, they have a 
positive real part of the dielectric constant). In this case, it is well know that 
no SPP resonances can exist. 

The results of this paper support those of Ref.4, although the analysis 
there, based on the angular spectrum of plane waves for the scattered field 
above and near apertures, attributes the enhancements to the contribution 
of a broad bundle of evanescent waves of this angular spectrum. As far 
as we know, that theory, which is also scalar, has not been tested with 
any exact procedure. We agree with the experiments and magnitude of the 
enhancements obtained in [4], as well as with the lack of contribution of SPPs 
to these peaks, as discussed there. We believe that our analytical calculations, 
validated by the FDTD simulations , help to support this agreement, even 
though it is not our purpose to enter in a direct comparison of our theory with 
that of [4]. We believe, nevertheless, that the mechanisms discussed here are 
sufficient to understand these enhancements as far as slits are concerned, in 
the same way as a recent publication done with the same procedure (cf. [6]) 
explains the enhancements for arrays of holes. In both cases: holes and slits, 
the characteristics of transmission of just one of these apertures determines 
the resultant properties of the array. However, whereas the properties of 
the transmitted modes in each slit for either s or p-polarization rules the 
enhancements of the array as discussed in Section 2; in the case of circular 
holes, the morphological resonances of each hole are a crucial factor that 
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contributes to the transmission peaks in addition to that of the Rayleigh 
anomalies of the arrays. 

We hope that the arguments presented in this work on the physics of this 
effect will shed light on discussions as regards the nature of these enhance- 
ments. 

This work has been supported by the DGICYT and European Union. We 
thank K. Ponizovskaya for assistance with the computer calculations. 
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Figure Captions 

• Figure 1: Scattering geometry. Period d = 3.5/zm, slit width a = 
0.5/im; the thickness of the slab h varies between 0.2/im and 4/im. 

• Figure 2: p-polarization. Transmittance versus illumination wave- 
length for the Ag grating parameters: period d = 3.5/im, slit width 
a = 0.5yum. Solid line: FDTD calculation. Broken line: Theory cal- 
culation based on Eqs.(4)-(6) (a) Thickness h = 3/im; (n e ff = 1.2 in 
theory calculation). (b) Thickness h = 2/im; {n e ff = 1.4 in theory cal- 
culation), (c) Thickness h = 1.2/xm; (n e // = 1.7 in theory calculation), 
(d) Thickness h = 0.6fim; (n e ff = 3.1 in theory calculation). 

• Figure 3: (a) p-polarization. Transmittance versus illumination wave- 
length for the Ag grating parameters: period d = 3.5/zm, slit width 
a = 0.5/zm, thickness h = 5fim. Calculation based on Eqs.(4)-(6) 
with n e ff = 1.2. (b) s-polarization. Transmittance versus illumina- 
tion wavelength for the Ag grating parameters: period d = 3.5/im, slit 
width a = 0.5/im, thickness h = 0.6/im. (c) p-polarization. FDTD 
calculation for the same grating as in Fig.2(d), but now the material is 
Cr. 
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